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Abstract 

Within a general theoretical framework we study the effective, deformation-induced interaction 
between two colloidal particles trapped at a fluid interface in the regime of small deformations. In 
many studies, this interaction has been computed with the ansatz that the actual interface config- 
uration for the pair is given by the linear superposition of the interface deformations around the 
single particles. Here we assess the validity of this approach and compute the leading term of the 
effective interaction for large interparticle separation beyond this so-called superposition approxi- 
mation. As an application, we consider the experimentally relevant case of interface deformations 
owing to the electrostatic field emanating from charged colloidal particles. In mechanical isolation, 
i.e., if the net force acting on the total system consisting of the particles plus the interface vanishes, 
the superposition approximation is actually invalid. The effective capillary interaction is governed 
by contributions beyond this approximation and turns out to be attractive. For sufficiently small 
surface charges on the colloids, such that linearization is strictly valid, and at asymptotically large 
separations, the effective interaction does not overcome the direct electrostatic repulsion between 
the colloidal particles. 

PACS numbers: 82.70.Dd; 68.03.Cd 

Keywords: Colloids; surface tension and related phenomena 
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I. INTRODUCTION 



The self-assembly of sub-/im colloidal particles at fluid (e.g., water/air or water/oil) 
interfaces has gained significant interest in view of various basic and applied issues such as 
;he study of two-dimensional melting investigations of mesoscale structure formation 
2I, and engineering of colloidal crystals on spherical surfaces jsl. The colloidal particles are 
trapped at the interface if the fluid phases wet the colloid only partially; this configuration is 
stable against thermal fluctuations and it appears to be even the global equilibrium state, in 
accordance with the experimental observation that the colloids immersed in the bulk phases 

n 

are attracted towards the interface 

In order to prevent coagulation, the colloidal particles are electrically charged. The ensu- 
ing repulsive force is well understood and at large separations d it varies like a dipole-dipole 
interaction oc l/c?^, because the monopoles vanish due to screening by counterions in water 

nn 

|5|. Nevertheless, several experimental findings have led to postulating an attractive ef- 
fective force between such particles with a range much larger than that of van-der-Waals 



forces 



, [lO|, lUl, I12I, [l3|, [ij, [l5| (but see also Ref. Il6l ). Spherical particles (radii 



R = 0.25 . . . 2.5 fim) at flat water-air interfaces exhibit the spontaneous formation of compli- 
cated metastable mesostructures consistent with the presence of a minimum in the effective 
intercolloidal potential at separations d/R ^ 3 .. .20 and with a depth of at least a few ksT. 

Until now, no unequivocal explanation for the appearance of these relatively long-ranged 
attractions is available. One possibility, which has been explored intensively in previous 
years, consists of an attraction mediated by the deformation of the interface (see Fig. [T]). 
This is similar to the so-called flotation force attracting particles floating at the surface of 
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But gravity plays no role 



water which is deformed by the weight of the particles 
for micrometer sized particles as described above. Instead, the electrostatic field around the 
charged colloids deforms the interface and gives rise to effective, capillary-induced interac- 
tions. There have been contradictory results about the properties of this effective interaction. 



In Refs. 



12 



19 



it has been argued in favor of an attractive 



brce decaying like 1/d (similar 



to the gravity-induced flotation force). In Refs. l20|, l21|, |22|, |23!, l2J this line of thought was 
shown to be invalid because the decay turns out to be much faster if the physical system 
consisting of the colloidal particles plus the interface is mechanically isolated, i.e., if the total 
force acting on this system vanishes in the limit of a macroscopically extended interface with 
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negligible borders. This is in principle the case for the experiments conducted in Langmuir 
troughs with lateral extensions several times the capillary length of the interface (i.e., orders 
of magnitude larger than the colloid radius). In this case the electrostatic force pulling on 



the interface is counterbalanced exactly 
water. In particular, the authors of Ref. 
However, this was in turn corrected in Re: 
totic decay cc 1/dJ was derived. In Ref. 122 



DV the electrostatic force pushing the particles into 



20 argue for an attractive force decaying like l/dJ . 



:s. 
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22 



and a repulsive force with an asymp- 



it was noticed that this latter result is actually 
unreliable because the linear superposition approximation employed in these calculations 
is not valid if mechanical isolation holds. Furthermore, since it is the electric field E, and 
not the electrostatic pressure oc E^, which obeys a superposition principle, the superposition 
approximation will be unsuitable if the main contribution to the electrostatic pressure stems 

with the conclusion that 



25 



26 



from cross terms in E^. This latter case was studied in Refs. 
the capillary-induced force is attractive and decays like l/c?^, which is the same asymptotic 
behavior as the direct dipole-dipole repulsion; whether the capillary attraction overcomes 
the electric repulsion must then be determined by a detailed analysis of the electrostatic 
problem. The total force is repulsive in the regime of small deformations of a flat inter- 
face (equivalent to the regime of small colloidal charges), for which the calculations were 



carried out 



25 



261]. Only for sufficiently large colloidal charges the capillary attraction 



may asymptotically overcome the electrostatic repulsion leading to a minimum in the to- 
tal effective potential 25|. However, further calculations — going beyond the linearization 
assumptions in treating the electrostatic pressure on the interface and the energy of the 
deformed interface — are necessary to substantiate this claim. The experiment described in 
Ref. [l^ is peculiar in the sense that the unperturbed interface is actually that of a relatively 



small spherical drop 
studied in Refs. [22 , 



et pending from a plate. T 



23 



27 



28 



29. In Refs. 



le importance of this finite-size effect was 



27 



29 it was found that the flotation-like 



decay l/rf is present because the plate breaks the condition of mechanical isolation, but it 
is quantitatively too small to explain the experimental observations^. 



Here we follow the approach of Ref. 



22 



in order to calculate the capillary-induced effective 



interaction beyond the superposition approximation in the case that there is a pressure fleld 



^ In this respect, Ref. 1 2 21 is incomplete but the same conclusion about the irrelevance of this effect is reached. 
Ref. |28| contains several important mathematical errors 29| and the corresponding conclusion erroneously 



disagrees with the actual irrelevance of the finite-size effect. 
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air / oil 




FIG. 1: Schematic drawing of the deformation of a fluid interface by electric fields due to charged 
colloidal particles trapped at the interface. Counterions gather on the side of the electrolytic phase 
(water) and a pressure field arises which pulls on the interface and on the particles. 



of general form 11 (r) acting on the interface in the limit of small deformations. The limit of 
small deformations corresponds to an analysis to leading order in the small dimensionless 
parameters ip and fn? defined in Eq. ([T]) below, which are measures of the force acting on 
the colloidal particles and the interface, respectively. In this manner, inter alia we provide a 



mathematically sound derivation of the results reported in Refs. |25, 



26l . In Sec. [TTl we derive 



the exact leading-order expressions for the capillary potential and show in which respect 
the superposition approximation becomes inconsistent if mechanical isolation holds. In such 
a case, the final result depends on the stage at which the approximation is introduced and 
applied. In Subsec. Ill Bl we derive the leading asymptotic dependence on d of the effective 
interaction using only some rather general assumptions on the form of n(r). The result 
for the effective potential energy is summarized in Eq. fl27|) . In Sec. IIIII we analyze the 
particular case that the pressure field n(r) is due to the electric field created by charged 
colloidal particles in a mechanically isolated system. We consider two limiting cases in order 
to solve the electrostatic problem: (i) water as a perfect conductor (i.e., vanishing Debye 
length), and (ii) the colloidal particles as point-like objects (i.e., vanishing particle radius). 
In both cases we recover the conclusion that the effective interaction is attractive but can 
overcome the electric repulsion only iiip > 1, which is outside the small-deformation regime 
considered here. Finally, in Sec. [IV]we summarize our results and discuss their relevance in 
connection with the experiments described in the literature. 
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II. FREE ENERGY OF EFFECTIVE CAPILLARY INTERACTION 



A. Exact results 

We consider identical^ spherical particles trapped at a fluid interface (see Fig. [2] for 
N = 2). We define the reference configuration as a flat interface in the plane z = and 
the colloids at a height such that the colloid-interface contact occurs at Young's contact 
angle 9 G (0, n). In this configuration, Sa denotes the circular disk delimited by the contact 
line on the colloid a, dSa is the corresponding contact line (of radius tq = Rsin6 for a 
particle of radius R) traced counterclockwise when viewed from the top'^, and S^en is the 
fluid interface (with surface tension 7), enclosed by a boundary Cl of typical size of the 
order of L representing, e.g., the vessel containing the system. The relative lateral positions 
of the colloids are kept fixed and thus we consider only vertical displacements ^(r = {x, y)) 
of the fluid interface and of the height Aha of the center of the colloid a relative to the 
plane z = 0. In the reference configuration there is a (vertical) force Fa acting on the colloid 
a and a (vertical) force per unit area n(r) on the meniscus. We define the dimensionless 
forces 

■■= -TT^, := / dA fl, (1) 



2vr7ro 27r7ro js^^^^ 

where dA is the element of interface area. The reference configuration is the equilibrium state 
in the absence of forces. Fa = 0, II = 0. Within the approximation of small deviations from 
the reference configuration, the free energy of the system with respect to this configuration 



is |22| 

:F = -f [dA Uvu\^--Ilu f de[Aha-u]^ -FaAha\+0{eF,en)\ (2) 

"'S'mcn - ^ . a=l ^ J dSa ) 

where di is the element of arclength. The free energy contains a contribution from the 
change of the contact area between the phases (two fluid phases and the solid particles), 
and a contribution from the work done by the forces Fa and fl via displacements from the 
reference configuration. 



^ In actual experiments there is a certain degree of polydispersity in size, shape, electric charge, etc., which 
we do not expect to alter the conclusions qualitatively if it remains small enough. 

Which fluid phase is to be the top one depends on the experimental setup. Usually gravity breaks the 
up-down symmetry. 
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FIG. 2: Top view (plane z = 0) of the reference configuration with two colloids, d is the fixed 
lateral distance between the colloid centers projected onto the plane z = 0. Si and S2 are disks of 
radius ro, the corresponding circumferences (counterclockwise) are dSi and dS2- The (projected) 
interface is ^men = 5R^\(5'i (J 5'2)- The position of any point on the plane is denoted with r; in 
particular, is the position of the colloid a. 



The values of u{r) and Ah^ in the equilibrium state are determined by minimizing this 
free energy. This leads to the following equations: 

Aha = {u)a - roEF^, ■= 7r~ f (')) (3a) 



27rro 



V'm = -in, re Smen (3b) 

7 

n„ ■ Vu(r) = Ef^ H , r G dSa, (3c) 

u{r) =0, re Cl, (3d) 

where is the unit vector in the outward normal direction of dSa- Equation fl3a|) is a 
geometrical relationship, Eq. (]3bp describes local mechanical equilibrium (the pressure 11 is 
compensated by the curvature-induced interfacial tension), Eq. fl3c|) describes mechanical 
equilibrium of the particle (the force Fa is balanced by the interfacial tension exerted at the 
contact line), and Eq. ( I3d[) represents a boundary condition at the external border Cl (for 
simplicity we take a pinned interface, but of course other physically reasonable boundary 
conditions are possible, the details of which are actually irrelevant in the limit L — oo 



we shall consider 



22l|). The free energy functional in Eq. ([2]) evaluated at the equilibrium 
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configuration can be simplified by using the relationships in Eq. ([3]) and we obtain two 
equivalent expressions: 

1 r ^ 

^eq = -^1 I dA \Vu\' + 7rjJ2[{u)l - {u')^-rlil^] (4a) 

If. ^ 
= "2 ldAIlu + 7r'j^roiF^[{u)a-roiF^]- (4b) 

From this point onward, we shall consider the particular case of two identical colloids. 
We assume that any external force acting on the system (e.g., gravity) is independent of 
the positions of the particles^, so that symmetry arguments will allow one to simplify the 
expressions. We use the notations ^^(r) := u{\r — rd). Hair) := n(|r — Tq,!), ep '■= ipaj 
Smcn,a '■= ^'^\Sa foT the Corresponding quantities in the presence of a single colloid located 
at position Tq,. This means that the function Ua{r) {a = 1,2) satisfies 

V Ua = Hq,, r G 5'men,a; (^^) 

7 

Ua ■ VUa{r) = Ef, T G OSa, (5b) 

u«(r) =0 reCL. (5c) 

In terms of these single-colloid solutions the configuration in the presence of two colloids 
can be written without loss of generality as 

n = Hi + + 2n^, (6a) 

U = Ui + U2 + Um, (6b) 

iF = eF + em, (6c) 
where we have introduced Bf '• — ^Fi — 

reflecting the symmetry of the problem. The 
fields Um{r) and nm(r) and the quantity Em introduced this way represent the corrections 
to the so-called superposition approximation, which is defined by setting 

Um = 0, (7a) 

= 0, (7b) 

£m = 0, (7c) 



^ This seems to be actually a reliable approximation for the experimental setups considered so far in the 
Uterature, but see the discussion in Sec. II VI 
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i.e., the effects of other particles on the single-particle configuration are neglected altogether. 
From Eqs. ([3], [5]) one can derive the following equations linking Um(r), nm(r), and Sm (with 
a,/3=l,2): 

V Um = r G 5'men; (S^) 

7 

n„ ■ - ^^^^^^^ = e„ - n„ ■ + v e dS^ (P^a), (8b) 

M™(r) = 0, r G Cl. (8c) 

The superposition approximation is violated even if 11^ = and Em = because of the 



boundary conditions at the contact lines 
non-vanishing 11^ was addressed in Refs. 



'Eq.jj8bl)). as pointed out in Ref. 



25 



26. 



22. The case of 



The capillary-induced effective interaction energy is defined as Knen('^) '■= -^eq — 2jF(,q, 
where 2J-'eq is the sum of the equilibrium free energies of the single-colloid configurations, 
i.e., for d oo. Vmcnid) depends parametrically on the (lateral) separation d of the colloid 
centers in the reference configuration. This is not the total interaction potential, which must 
include, e.g., the direct electrostatic repulsion between charged colloids, not considered in 
the expression ([2]) for the free energy. The effective interaction energy can be written as 
Knen = Kup + Korr, whcrc V^np is the rcsult of imposiug the superposition approximation 
(Eq. (I?])), and Korr is the correction to this approximation. If we use the expression ( llal) . 
we obtain 

"K.p = —f [dA (Vixi) ■ (Vm2) +-f fdA I VmsI' - 27r7 {uj - {u2)l), , (9a) 

"V^corr = -\^ L ^^^rn) " V(M™ + 

^ 5'mcn 

- 27r7 {UmiUm + 2^2) - {Um)l{Um + 2M2)i)i " 27r7r^£:„ (25^7' + Em)- (9b) 

On the other hand, using expression ( ]4b[) . we obtain 

Vsup = - dA U1U2 + dA U2U2 + 2njrQEF{u2)i, (10a) 

5'mcn ^ Si 

^Vcorr = - dA [Ea Um + 211^ U2 + Um Um] 

^ Sm.cn 

+ 27r7ro((£F + £m)Um + £m{Ui + U2))l - 27ljrlEm{'2EF + ^m)- (lOb) 
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We emphasize that the two alternative expressions of Vsnp or V^orr are not equivalent (but 
their sum Vmen is), and they in turn differ from V^up as computed in Ref. (see Eq. (39) 
therein), which was derived by inserting the superposition ansatz directly into Eq. ([2]) here. 
Application of Gauss' theorem with Eq. (jSaj) leads to 

'^Hup - 'V^up = 2vr7 (roiu^ + «2)|^ - (^2 - (^2)1)'^ = %orr - "Korr. (H) 

The superposition approximation is inconsistent asymptotically in cases in which "Vsup— *K)up 
doesjiot decay more rapidly than '^V'sup as function of the separation d. As remarked in 



Ref. 



22 



there are indeed relevant cases in which this consistency condition is not fulfilled 



[see Subsec. IIIBp . 



B. Effective potential in the intermediate asymptotic regime tq <ti d <^ L 

In this subsection we compute Knen(c?) asymptotically in the intermediate range tq <C 
d <^ L. For this purpose, we have to make some restricting assumptions which, however, 
seem to be satisfied in the experimental setups investigated so far. First, in view of the 
discussion in the Introduction concerning the electrical fields, we assume the proportionality 

C-iniHsl. (12) 

This is valid if the interface stress is quadratic in a field satisfying linear superposition in the 
two-particle configuration. Examples are given below by some specific electrostatic models 
(see, c.f., Eqs. ( l3Ti) and ( H9l) ). Second, we assume that the single-colloid pressure 11 decays 
far from the colloid as 

n(r) ~ r"", n>A. (13) 
In the experimentally relevant case of charged particles at water interface, it has been estab- 



lished both theoretically {4] and experimentally 30j that n = Q. (This can be understood 
easily: the charge of a particle induces a screening image charge in the water, so that the 
distant electric field is dipolar.) The constraint n > 4 will allow us to estimate the integrals 
appearing in Eqs. [TUI) by approximating them by the contribution of the regions near the 
colloids. This condition excludes, however, the effect of an external electric field (in that 
case the resulting pressure does not have to decay at all) and the case that the colloidal 
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charge is not perfectly screened so that the distant electric field corresponds to a monopole, 
i.e., n = 4 (this would occur if both fiuid phases are dielectric, e.g., air and insulating oil). 

The quantity iu — '2iF is the net (vertical) force by an external agent acting on the total 
system consisting of two colloids plus the interface (see Appendix This can be, e.g., 
gravity (if the colloid is large enough for it to be quantitatively relevant), dispersion forces 
by a substrate closely beneath the interface (this effect can be modeled similarly as gravity, 
see Appendix |B]), or an optical tweezer pushing the colloid vertically. Since we have assumed 
previously that this external force is independent of the positions of the particles, it is given 
by the sum of the net forces in the single-particle configuration (i.e., if they are infinitely 
far apart from each other): 

in-2iF = 2{en-ep). (14) 
From the definitions in Eqs. ([H [6al) we can write 

in = 2{en + Su,^ - eu) (15a) 

with 

:= / dA Um, := / dA U2. (15b) 

2T^iro Js^^^ 27r7ro 



For (i — > 00, it is clear that 



Si 



eu ~ ^. (16) 



In this limit, we note that receives its main contribution from the regions around 5*^ 
implying 

provided n > 4, because |nm(r)| ~ ^/I^dyuJJ^r^^^V^ in those regions which provide the 
dominant contribution to the integral. Therefore, from Eqs. (!6cl [T4l) one obtains asymptot- 
ically 

1 

em = - £12 ~ ~ (18) 

With these simplifying assumptions, we shall compute analytically the behavior 
of Vaien{d — > oo) to leading order in 1/d. More precisely, on dimensional grounds^ 



^ Due to the rapid decay of 11 far from the particles, in the limit L — > 00 the length scale L enters into the 
problem only if en 7^ ef and in that case just as an upper bound on d to avoid a logarithmic divergence 



(see the discussion after Eq. (jAl[) ). 
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the expansion parameter is ro/c?. In principle, fl can contain and thus introduce ad- 
ditional length scales, for example the Debye length if the fluid phase is an electrolyte. 
This complicates the problem, which then has to be analyzed numerically (see Subsec. [Ill Bp . 



The superposition approximation: Within the superposition approximation Knen('^) 



was computed in detail in Ref. 



221 . Here we sketch briefly the estimate of the asymptotic 



behavior of ''V^up (compare the three terms in Eq. fllOal) ): 

2/" dAu2-^7irl€u + T^rl{u2)2 

^ '5'mcn.2 

because the main contribution stems from the regions around Sa-, and 



(19a) 



f dA XI2U2 ~ n(c/)M(c/) / dA 

5i isi 



after expanding around r2 = d. Thus 



uid) - j^rlU{d) 



2 1 dAu2 + Trrl{u2)2 



(19b) 
(19c) 

(19d) 



because Il{d) is asymptotically subdominant compared with u{d) (see Eqs. (IT3 | lAip ). The 
difference "Vsup — ^ V^up in Eq- f|TT1) between the two implementations of the superposition 



approximation can be estimated as 



22I (see Eq. (IXB ^ 



(^^1 + ^2)^y ~ -ro(Mi)iV\(d) oc n(rf) ~ — , 



{{u2-{u2Wf),^-r',\Vu{d)\^^ 



(20a) 
(20b) 



l/c/2(«-i), iieF-eu = 0, 

Thus, whenever ep — £u 7^ (i.e., the system is not mechanically isolated in the sense 
that there must be a force acting on the boundary Cj;, of the interface to compensate this 
non-vanishing net force) one finds 



'Vsnp{d) ~ 27r7ro(eF - £n)u{d) ~ 27r7ro(eF - £nf In y + const 



(21) 



corresponding to an attractive force irrespective of the precise form of the function n(r). 
(The additive constant, which does not affect the physical conclusions, depends on the 
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precise form of the boundary condition at Cl.) Physically, ''V^up is the work done by the net 
force 27i'y{eF — upon a vertical shift of the subsystem consisting of one colloid plus its 
surrounding interface (behaving like an "effective particle") by an amount u{d) due to the 
deformation induced by the second colloid. In this case (i.e., ep — £u 7^ 0), the difference 
"'Vsup — ^Vsup decays more rapidly than ^V^up('^) and both expressions "Kup and ^Vs^p agree 
asymptotically. Equation (!2T!) exhibits the same dependence on d as the potential energy 
associated with the flotation force; this is discussed briefly in Appendix [Bl 

If ^i? — en = (corresponding to mechanical isolation), '^Vsup{d — » oo) ~ 1/d"- (see 
Eq. (1131) ). This decay agrees with previous findings I21L I22I ] , but the reliability of this 



result is unclear because the difference "Kjup ~ ''Vsup decays with the same power law. As 
a matter of fact, the amplitudes of the asymptotic decay of "Kup and ''l^up differ and are 
in turn different from Eq. (52) in Ref. |2J, because there actually neither of the two repre- 
sentations of JFpq was used. However, for the leading behavior this is unimportant because 
Vmen is asymptotically dominated by the correction to the superposition approximation. 

Beyond the superposition approximation: If the interface deformation field in Eq. (I6bll 

is evaluated near colloid 1, the term U2 + Um is dominated by U2 in the absence of mechanical 
isolation and by Um in the case of mechanical isolation (see, c.f., Eqs. (lAlt IA9P ): 

[ Uynid) ~ c/~"/^ if eF-eu = 0, 

U2+Um ^ ( (22) 

I U2{d)^\nd, iiep — EuT^O- 

Thus the superposition approximation holds if and only if the system is not mechanically 
isolated. Otherwise, the correction Um is dominant and, for asymptotically large separations 
d, Vuien ~ VcoTT- Hence in the following we take ep — £u = 0- Since the main contributions 
to the integrals stem from the regions around Sa one obtains the estimates 

- [dA [U2U^ + 2UmU2 + IlmUm] / rfA [H2 + SH^ Mg] ~ . (23) 

Using the asymptotic decay of Um given by Eq. (1A9P we find 



b 



Korr(c?) 



- I dA [H2 Um + 2H^ U2] + 2n-fro {epUm + £11^^1)1- ^n-frlepenm ~ 

(24) 

because each term scales like c?""/^ and there is no reason for mutual cancellations. In this 
expression one can identify two distinct contributions with a simple physical meaning. After 
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inserting Eqs. (l6bl [6cl) into the relationship (l3al) we write A/i2 = A/i2 + A/i^ with 

A/i2 := (Mi+M2)2-roeF, (~ (^2)2 - ro + (25a) 
A/i^:=K)2-ro£^ (25b) 

where the asymptotic decays are giving by Eqs. (JTS llAQl) . Accordingly, expression ( l24l) shows 
that ^Korr is dominated asymptotically by^ (i) the work done by the additional pressure 
2Ilm if the single-colloid configuration is deformed relative to the reference configuration, 
that is, — / dA 2IlmU2 + 27['~frQe-[i^Ah2, and (ii) the work done by the forces acting in the 
single-colloid configuration upon the additional displacement Um, that is, — J dA 112 Mm + 
2']TjroeFAhm. 

The asymptotic behavior of ""Vcorrid —>■ 00) can be derived by estimating the behavior of 
the terms in Eq. ( ]9b[) individually, as in the case of Eq. (flObll . with the result 

"Korr(rf) ~ -27 f dA (Vu„0 ■ (VU2) - ~ (26) 

The asymptotic decay predicted by Eqs. ( 12^ and ( l26l) must agree both with respect to the 
decay exponents and the amplitudes, because the difference "Vsnp Kup decays asymptoti- 
cally more rapidly (see Eq. (!20!) for — en = 0). The sign of the force described by Vcorrid) 
is not evident from the outset, but in the applications we shall consider later, it turns out 
to be always attractive. 

It is interesting to compare our result with the corresponding one in Ref. |26). After 
noting the equivalence Um ^ 2ui2 in the notations, one finds that Eq. (7) in Ref. l26| is 



identical with the integral term in Eq. (!26l) here. We obtain an additional term oc epEum 
because we treat the contribution by the colloid to the free energy functional in full detail, 
while in Ref. 2<^ the point-particle approximation (ro — * 0) is used from the outset (compare 
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our additional term 



Eq. fHa|) here with Eq. (1) in Ref. '26^. As a consequence, in Ref. 
is lost as a singularity of the integral term, which is regularized there by introducing an 
unknown cutoff length expected to be of the order of ro (see also Refs. l2l|, 122). The analysis 
of realistic models in Sec. IIIII will show that the quantitative contribution of the term 
oc epSUm to effective interaction V^en is actually larger than but proportional to the 
contribution from the other term in Eq. (!26l) . 



Note that by symmetry (^2)2 — (1*1)1 and (wm)2 — {um)i 
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To summarize, the capillary-induced effective interaction between two colloids is given 
by (compare Eq. f[T3|) defining n) 

Knen(rf) ~ < « d) . (27) 

[ Vsnp{d) ~ Ind, if Ef -Bu^ 0, 

III. APPLICATIONS 

In this section we compute Vmen('^) for e^? = en (mechanical isolation) for different re- 
alistic models of 11 derived from the solution of the electrostatic problem within various 
approximations. We note that the asymptotic decay of Vjncn{d) is the same as that of the 
direct electrostatic repulsion, so that these detailed calculations beyond the asymptotic anal- 
ysis of Subsec. Ill Bl are necessary in order to be able to address this fine-tuning problem and 
to determine whether the total force is asymptotically attractive or repulsive. 



A. Ideally conducting fluid phase 

The simplest model consists of approximating water by an ideal conductor. Formally 
this corresponds to the limit of zero temperature, so that the osmotic pressure of the mobile 
charges accumulated at the interface vanishes and the Debye length is zero (see, c.f., 
Eq. fHTj) ). In this case the electric field E = Ee^ is always normal to the interface and the 
pressure is given by Maxwell's stress tensor evaluated at the insulating side of the interface 
(we use Gaussian units), 



n = — e, • ( eEE - -eEH] ■ 
47r V 2 



= ^E\z = 0^), (28) 
=0+ 8^ 



where ei is the dielectric constant of the insulating phase [z > 0). 

In the present context the electrostatic problem of a charged sphere partially immersed 
in a conducting fluid has been solved numerically and semi-analytically in Ref. |3l|. There it 
has been shown that the single-colloid pressure field exhibits an integrable divergence upon 
approaching the three-phase contact line and that asymptotically it displays the familiar 



dipole behavior. The following approximate parametrization (Eq. (1.4) in Ref. 
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expressed 



in terms of our notation) incorporates these properties and is sufficiently accurate for our 
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present purposes: 



n(r) = ^ bif^) 











r 


- 1 




r 


Jo 









6(/i):=-/i(/i + l)(/i + 2)(/i + 3), (29) 



where < /i < 1 is a fitting parameter the precise value of which depends on the contact 
angle 6 and the dielectric constant ei. This expression is normalized so that en = £f and 
it corresponds to an exponent n = 6 independent of the choice for /i. Therefore, far from a 
colloid, the single-colloid electric field |E| = ^/SnTl/ei is that of a dipole perpendicular to 
the interface, the strength of which is given by 



\p\= hm e,r'\E{r)\ = e, rg J ^^^^^^^^^ (30) 

In the presence of two colloids, we take E ^ Ei + E2. This approximation allows a 
simplification of the calculations and should not alter the physical picture significantly. The 
approximative character is due to possible violations of the electrostatic boundary conditions 
at the surfaces of the colloidal particles: The additional polarization of colloid 1 induced by 
E2 will actually lead to an electric field E = Ei + E2 + 5E in the neighborhood of colloid 
1, with an induced electric field |5E(r, (i)| ~ x(i")-^2('^), where the "effective susceptibility" 
x(r) is expected to be at most of order unity. Thus, our conclusions will be qualitatively 
correct with a quantitative error of a factor of order unity. Under these conditions the field 
nm(r) defined by Eq. fIBal) is given within this approximation by 

n^(r) = Vn(|r-ri|)n(|r-r2|). (31) 

The integrals in Eq. fl26|) are computed by using the expressions in Eqs. (I A 11 IA9p so that in 
the limit vq d one obtains (for details see Appendix [U]) 

- r-^y (32a) 



11+ 1 \ d / 

and 



[ dA {VuJ ■ (Vu,) ^ ^^M(/i) rlel (^) ' , (32b) 



with 

M(/i) := g (/i + !)&(/") J^dv 2^1(1 - /i, 4; 5; v) 2Fi{^^, 1; 2; v) (33) 
and 2-P'i («,/?; 7; -z) is the hypergeometric function (see Subsec. 9.1 in Ref. 132)- We have 



checked numerically that M(/i) ^ /i/5 with a maximum error of ~ 1.3% within the range 
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< /i < 1. The final result reads 

V^M ^ ""V^orM ^ -^^^[1 + Mif,)] ^rl 4 (^)' . (34) 

This corresponds to an attractive force. We note that the dominant contribution to Vmenid) 
stems from the term proportional to in Eq. (l26l) . the more so the smaller /i is. Smaller 
values of fi correspond to an increasing importance of the electric field near the colloid (see 

Eq. (EHD). 

Vuien{d) is to be compared with the potential energy due to the direct electrostatic repul- 
sion Vj-ep{d) of the charged colloids. The potential energy of one dipole in the field of another 
identical dipole is p^/(eic?^) so that Vrep{d) ~ /{2eid^) for large d. (One must divide by a 
factor 2 because within our model no work is done on the image charge inside the conducting 
phase forming the dipole.) Collecting the results, we find that the total interaction energy 
at large separations is given by 

/r„N 3 r A 1 

(35) 



Vtotal(c^) = Vrcp{d) + Vmen(rf) ~ ATX^irl E fK^) ' 



d 



1-^^[1 + M{^j)]ef 



Hence we see that the attractive capillary potential is proportional to and the direct 
electrostatic repulsion is proportional to Ep- Thus in the limit Ep 1, on which our calcu- 
lations are based, the electrostatic repulsion is always larger than the capillary attraction. 
The leading-order analysis of this model predicts an attraction only if the charge of the 
colloid is large enough so that Ep = 0{1). The critical value above which there is attraction 
is given by 

'^'-^^^^ = 4[1 + MM] ' ^^^^ 
which lies in the range 1/4 < Ep^crit < 5/12 for < /i < 1. 



B. Finite Debye length 

Consider now the more general case of an upper insulating phase (dielectric constant ei) 
and a lower electrolytic phase (dielectric constant €2, electrolyte concentration wq) at a finite 
temperature T. The stress tensor acting on the interface is due to the difference of Maxwell's 
stress tensor just above and below the interface, plus an osmotic pressure posm by the excess 
of ions concentrated close to the interface: 

n = — ■ ( eEE - -eE^ 
An \ 2 



z=0+ 

+ Posm- (37) 



z=0- 
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With ^(r) denoting the electrostatic potential at the interface, Ez^±{r) := ■ E(r, z = 0^) 
as the normal component of the electric field at the interface (with eiEz^+ = e2-Ez,_) and 
E|j(r, z = 0) as the (continuous) parallel component at the interface we have 

EE - -eEH ] ■ e. 



An 



£2 - ei 



2=0- 



61T 



ei ^2 



-E: 



£2 



E 



(38) 



and with An(r) as the excess ion concentration at the interface 



Po 



kBTAn = noksT 



exp 



+ exp 



6$ 



(39) 



assuming that the electrolyte is dilute and consists of monovalent ions. 

In order to solve the electrostatic problem we introduce two simplifications: We apply 
the Debye-Hiickel approximation for dilute electrolytes and we approximate the extended 
colloid by a point charge q at its center, i.e., we retain only the monopole term of the source 
of the field. This 'monopolar' approximation corresponds formally to the limit of vanishing 
contact radius, tq 0, and it can be expected to provide the correct field at distances from 
the colloid which are large compared to tq. This is in fact complementary to the limiting 
case KTo ^ 1 considered in the previous subsection. From Eq. ( |39l) we find 



Po 



where the screening length is given by 



in 



(40) 



e2kBT 



(41) 



For pure water €2 = 81 and Uq ^ 10~^ mol/1 (that is, pH = 7) lead to a screening length 
^ 1 yum at room temperature. 

We consider first the case of a single colloid corresponding to a point charge q located 
at the flat interface between an insulator and an electrolyte. Proceeding along the lines of, 



e.g., Refs 



m 



we find: 



E. 



$(r) = ^J„(«:r), 
e2r 



E|i(r) 



e2r 
2q 



dr 



t2r 



2 2^c('f") 



(42a) 
(42b) 
(42c) 
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where is the unit radial vector pointing away from the center of the colloid, and the 
auxiliary functions X„(/c) are given by integrals over a Bessel function: 



X 



dx Jo(x) — , , 

+ 



1:2 



— / dx Joixj , 



Xa(A:) 



(43a) 
(43b) 
(43c) 



These integrals have been computed numerically; the analytic expressions for their asymp- 
totic behaviors are derived in Appendix [Dl Thus, taking e2/ei ^ 1 (the ratio e2/ei is 
approximately 81 for water in contact with air) in the intermediate asjTiiptotic regime 
1 <^ Kr <^ (£2/^1)^ one obtains the following expressions for the potential and the field 
components (see Eq. fH2l) ): 

1 



$(r) 
Eiifr) 



2q_ 
2q_ 

2q_ 

€2^2 



_e2 (nry 
1 

+ 

Kr 



+ e" 



1 !i («:r)3/2e— 



(44a) 
(44b) 
(44c) 



62 (Kr)2 

Asymptotically $(r) decays as 1/r^ and the electrostatic interaction energy q^{d) of a second 
charge g at a distance d from the first one decays likewise as 1/d^. This is the celebrated 
dipole repulsion between charged colloidal particles first conjectured in Ref. Q. However, due 
to the large value of the ratio ^2/^1, closer to the charge there is a crossover to a screened 
Coulomb potential $(r) oc exp (— /tr)/r. One can introduce a crossover length Tcross defined 
from the asymptotic behavior given by Eq. (jHj) as 



ei 



£2 ('^r cross 



(45) 



This equation has solutions only if ei/e2 < (S/e)^ ^ 0.54, in which case the relevant solution 
'^'"cross is larger than 2 and depends only weakly, i.e., logarithmically on the precise value of 
the ratio ti/t2- For t2/ti = 81 this gives /tr cross ~ 8-7. The parallel component Ey exhibits 
the same crossover, but the normal component always decays algebraically and is much 
larger than the parallel component. Nevertheless, the contribution of Ez^+ to the pressure 
field n(r) is reduced by a factor ei/e2 (see Eq. (I38l)). so that n(r) will also exhibit the 
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crossover at a distance r ^ rcmss- From Eqs. (138| [39]) one finds 



n(r) = 'jroK^ 



Bp 



where (see Eq. ([1])) 



£1 

£2 



£2 



£f — £n 



1 - - 1 J^(/€r) + ( 1 - - 1 lliKr) + (/€r)X'(«:r 



62 



27re27ro 



with the dimensionless function 



^00 

P(Kro) := / dx — 

J Kro 



£1 
62 



+ 



1 - I ) X,^(a:) 



-2 Kro 



£2 



(46) 



(47) 



(48) 



where the last line is the leading (zeroth order) contribution in an expansion in terms of the 
small parameter ei/e2- Equation ( H71) allows one to express the parameter q, the value of 
which is often uncertain, in terms of the more convenient parameter ep. Figure [3] shows how 
n(r) attains the dipole limiting behavior oc due to the electrostatic pressure oc E'^_^ 
beyond the crossover distance; at shorter distances, n(r) is dominated by the term oc Ey 
and, to a lesser extent, by the osmotic pressure. The figure indicates also that the main 
contribution to the total force eu stems from the regions close to the particle, as evidenced 
also by the formal divergence of V^^tq 0) (see Eq. (HHl) ). Thus, the precise value of en will 
be affected by the fact that our solution of the electrostatic problem within the monopolar 
approximation is expected to be reliable in principle only sufficiently far from the particle. 

We consider now two identical point charges at the interface separated by a distance d. 
Within the Debye-Hiickel and point-charge approximations, the solution of this electrostatic 
problem is given by the superposition of the single-colloid fields. In this case the field Tlrni^) 
in Eq. (l6a|) reads (as in Subsec. Ill At the subscript a{= 1, 2) denotes that the corresponding 
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K r 

FIG. 3: The dimensionless stress n := Ii/['^rQK?'ep /V^nro)] due to a single point charge at the 
interface for €2/61 = 81, as given by Eq. ()46p (thick soUd Une). The thin hnes represent each of the 
three additive contributions in Eq. (jl6]) : due to E^,, only given by X^; due to -Ey, only given by Xc\ 
and due to the osmotic pressure, only given by Xq. As expected, one observes a crossover in 11 at 
a distance r comparable with rcross defined by Eq. ([^5]) . 




field is evaluated with respect to particle a; see Fig. |2]for the notation): 



n^(r) 



£2 - ei 



57r 



ei 

€2 



(E,,+)l(E,,+)2 + (E|i)i-(E||)2 



57r 



27re2 K^(|r - ri||r - r2\Y 



— ( 1 - — ) Xb{K\v - ri|)J5(/«|r - r2|) + 

€2 V ^2, 



£2 



1 - — I XJ/tlr - ri|)Je(K|r - ral);^^ — ^ • — ^ + 



|r — ri| |r — r2| 



k^\y - ri||r - r2| Ja(«|r - ri|)2ra(«;|r - ra]) 



Therefore one has (see Eqs. (]15bl S 
where 



Vm{i^ro, Kd) 



Ef-i 



Vm{Kro, Kd) 



1 



(49) 



(50) 



(^1 - 1^ Jfe(xi)Jfe(x2) + (51) 



1 - — ) Ic{Xi)Ic{x2) "^^ + XiX2Xa(a;i)2a(a;2) 

£2/ a;iX2 
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in terms of x := nr and Ta). Dependences on tq and d enter through the 

specification of the integration domain S^aen- Using this expression and Eq. fH3|) . we have 
integrated numerically Eq. ([H]) using bipolar coordinates in a domain corresponding to a size 
L ^ 100/k. Taking this numerical solution and the exact expression in Eq. flAip . we have 
subsequently computed numerically the effective interaction energy as Knen('^) ~ ^Korrlc^) 
by using Eq. (llObp . As in the case studied in Subsec. IIII Al (see the text after Eq. (l34ll ). the 
last term in Eq. fllObI) turns out to be larger in magnitude than the remaining terms but 
they are approximately proportional to each other, so that for the purpose of understanding 
the numerical results we can consider the following proportionality (upon the application of 
Eq. 



Figure H] shows indeed that for two choices of ktq this approximation is reasonable at large 
d. Numerically we find that the term oc epSnm contributes ca. 70% of the total meniscus- 
induced potential V^cn- Asymptotically for large d it exhibits the predicted l/ci'^-behavior 
(Eq. (!27|) with n = 6). More precisely, the asymptotic behavior of Knen can be obtained 
as an expansion in terms of 1/d by utilizing the two-peak structure of Ilm (as for the 
general discussion in Subsec. IIIBp : One assumes that for d ^ 00 the main contribution 
to the integral in Eq. fl^T]) stems from the regions near the colloidal particles. The leading 
asymptotic behavior of Vm is given simply by the lowest-order term of a Taylor-expansion 
about Xi ~ ktq, X2 ~ i^d and Xi ~ nd^ X2 ~ i^tq in the integral in Eq. (15T|) with the 
asymptotic decay for the X-functions obtained in Appendix [DF: 



where 



Pm{i^ro) := 2 / dx 



00 



62 {nd)^ 



£2/ X 



(54) 



^ Next-to-leading terms in the expansion require a more elaborate calculation than the simple Taylor- 
expansion, which yields undefined expressions for them. 
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FIG. 4: The dimensionless capillary-induced potential energy Fmen := 10^T4nen/[7'^o^|'/'P('«?'o)] 
for two different values of hvq and the choice €2/^1 = 81. Thick lines correspond to the capillary- 
induced potential given by Eq. (jlObp (see Eq. (j27p ). whereas thin lines show the approximation 
oc ep^Um given by Eq. ([52]) . If kvq is large enough, a minimum appears. 



and the second line is the leading (zeroth order) contribution in an expansion in terms 
of the small parameter ei/e2, obtained from applying Eqs. (]D2l ID13I) . According to this 
expression, Pm{i^ro 0) = 4 and p(/tro) decreases monotonously with an asymptotic decay 
p{nrQ — > 00) ~ 2/ {hvq). This result for Vm{i^fo, f^d) means that the asymptotic l/(i'^-decay is 
determined by the normal component Ez{d) of the field and by the potential ^{d) stemming 
from the osmotic pressure. There is no contribution from Ey to leading order in l/rf due 
to the geometrical factor xi ■ X2 appearing in Eq. (!5T]) . Finally, from Fig. IHwe infer that if 
kvq is large enough, the asymptotic decay breaks down as a reliable approximation and a 
minimum appears in Knenl*^) at a separation d which is a few times k~^. 
The energy due to the direct repulsion of the two colloids is given by Eq. 

lai^d) 



{Kd){Kro)V{KrQ 

and the total energy is Vtotid) = Vrep{d) + Vmen{d). With the approximation (l52l) one has 



(55) 



V^tot(rf) 



KroP(Kro) 



Kd 



ep KToVmiK-ro, Kd) 



Asymptotically for Kd ^ 1 this expression reduces to (see Eqs. (153| IDS 



ei/e2 



ktqY V{Kro) V d 



^2, 



1 + — (K(i)^e '"^ - KroPmiKTo) ep 



(56) 



(57) 
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to be compared with the potential obtained in the previous subsection (Eq. (135|) ). As in 
that case, the capillary-induced potential Knen is reduced by a factor <^ 1 with respect 
to K-ep, and the total potential can be asymptotically attractive only if ep is above a critical 
value 

^F.crit = 7 T, (58) 

which to leading order is independent of the small ratio ei/e2 and turns out to be bounded 
as £:_F,crit > 0.38 (see Fig. [5]). Moreover, we note that even Hep ~ 1 the term in brackets 
in Eq. fl571) changes sign at a separation d comparable with the crossover length rcross (see 
Eq. This is a consequence of the crossover in Vj-epid) oc ^{d) and suggests that although 

the potential is asymptotically attractive, it reaches a minimum at a distance d ~ Tcross and 
turns repulsive for closer separations. This effect can only be enhanced by the deviation 
from the asymptotic 1/d^-decay in V^en('^) observed in Fig. H] for ktq ~ 1. 

These conclusions based on the approximation (l52l) are supported by the corresponding 
full numerical calculations. More precisely, from a fit to these numerical results in the range 
0.1 < KTo < 2.0 for the ratio €2/^1 = 81 we find that the critical value is given approximately 
by 

£F,crit(«ro) ^ 0.34 J 1 + . (59) 



Figure E] shows a plot of Vtot{d) for a typical value ep = 0.6, exhibiting a shallow minimum 
if it is asymptotically attractive. In summary, the capillary-induced attraction can domi- 
nate asymptotically only if ktq is sufficiently large so that ep > £F,cTit, in which case one 
necessarily has ep ^ 



IV. SUMMARY AND DISCUSSION 



We have carried out a thorough analysis of the effective interaction potential Knen('^) 
between two colloidal particles floating a distance d apart from each other at a fluid in- 
terface due to the deformation of the interface caused by the particles. The main result 
is summarized in Eq. (1271) for the asymptotic behavior of Knen('^) as c/ ^ 00. One flnds 
two qualitatively different cases depending on the (dimensionless) force ep acting on the 
particle, and the one, en, acting on the interface (see Eqs. [HI)): (a) if the system is 
not mechanically isolated {ep — eu ^ 0), the superposition approximation is valid and the 



24 



2 
1.5 
H 1 
0.5 


0.5 1 1.5 2 

Kro 

FIG. 5: The solid line gives the critical value e_F,crit as a function of Kro specified by the fit to the 
numerical results (Eq. (j59]) ). The dashed line corresponds to the approximate expression provided 
by Eq. (pi). 
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FIG. 6: The dimensionless total potential energy t4ot := 10^Vtot/[7?'o^F/'P(K?'o)] for different 
values of hvq and the choices 62/^1 = 81 and £f = 0.6. Note that the total energy is about 10 
times smaller than Vmen (the normalization of Vtot here differs from that of I^nen in Fig. H] by a 
factor 10 /ep) and exhibits a shallow minimum provided that Kro is not too small. 

asymptotic dependence Vmen{d 00) is universal in the sense of being independent of the 
detailed distribution of the forces which deform the interface. The physical reason for this is 
that each particle together with its surrounding interface can be considered as an effective 
particle, (b) In the opposite case of mechanical isolation {ep — £n = 0), the corrections to 
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the superposition ansatz are dominant and the asymptotic decay of Vmcn(c?) depends on the 
asymptotic properties of the pressure field acting on the interface. We identified two sources 
of violation of the superposition approximation: the boundary conditions of the deforma- 
tion at the particle-interface contact line are violated, and the pressure field acting on the 
interface does not satisfy a superposition principle. 

These results hold under the rather general assumptions that (i) the interfacial defor- 
mation is small {eF,£n ^ 1), (ii) the pressure field 11 = IIi + 112 + Slim exerted by the 
particles on the interface satisfies the scahng relations in Eqs. ( IT2| [T3!) . and (iii) the exter- 
nal force acting on the particles is additive (Eq. ( |T^ ). Assumptions (i) and (iii) are quite 
general and there is no evidence that experimentally they are not fulfilled. Assumption (ii) 
can be viewed as a consequence of the condition that the pressure II is derivable from the 
stress tensor of a Lagrangian quadratic in an underlying field, as it is typically the case for 
physical systems: Besides stresses induced by electric fields as addressed here, there can be, 
e.g., elastic stresses which arise if one of the fluids is in a nematic phase 33||. (This latter 
case provides also an example in which the asymptotic decay (!27j) is actually modified by 
geometrical constraints which complement the scaling relation ( fT2l) .) 

As an application for the case of mechanical isolation, we have considered the paradig- 
matic system of electrically charged colloids (see Figs. [1] and [2]), which is also of direct 
experimental importance. The pressure deforming the interface is due to the electro- 
static field emanating from the particles and gives rise to an effective attractive interaction 
Vmen{d oo) ~ 1 / . Siucc the direct electric repulsion between the colloids also decays 
\4ep ~ ^/d"^, this asymptotic analysis is insufficient to determine whether the total potential 
Knen + Kep dcscribcs an asymptotically attractive force as seemingly observed experimen- 
tally. To this end, as a model for the experiments it is necessary to consider in detail the 
challenging electrostatic problem posed by two charged spheres partially immersed in an 
electrolyte. Within this approach we have studied the electrostatic problem in two simplify- 
ing limiting regimes: (A) water as one of the fluid phases is a perfect conductor, and (B) the 
colloidal particles are replaced by point charges (monopolar approximation). Within both 
approximations we concluded that V^cn oc e^,, while V^ep (x ep, so that the direct repulsion 
dominates asymptotically; the capillary-induced attraction is predicted to dominate only if 
Ep is larger than a threshold value £:_F,crit ~ 1 (see Eqs. fl36 | [59]) ). i.e., the deformation of the 
interface has to be large, which is outside the range of validity of the present analysis based 
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on small deformations {ep ^ !)• 

Assumption (A) is a simple model of the experimentally relevant case of polystyrene 
or glass colloids floating at the interface between salty water and air or oil such that the 
Debye length of water is much smaller than the radius of the particle-interface contact 
line tq. Only the residual charges at the interface between the colloid and the air or oil 
phase determine ep. The charge density cr at this interface is rather small compared with 
the nominal surface density of dissociable groups; we note that in the case of polystyrene 
colloids the precise microscopic origin of this charge density is still unknown. By dimensional 
analysis, ep = (o"^ ro/ei^) G{ec/^i, 0), where G is a dimensionless function depending on the 
contact angle 6 and the ratio of dielectric constants of the colloidal particle and the insulating 
fluid, ec and ei, respectively. The electrostatic solution |31j] yields G{1,9) ^ 3/ sin'^ 6* for 



hydrophobic colloids. In Ref. |34r the single-colloid meniscus deformation was measured 
around rather large glass spheres {tq ^ 200 /im) at a water-oil interface (7 ~ 0.05 N/m, 
ei ~ 2), which are slightly hydrophobic {6 ~ 120°). Using Eq. ( lAlj) a value ep ^ 0.4 
was inferred from the measurements, corresponding to a charge density of a ~ 70yuC/m^ 
(or 5 X 10~^ e/nm^) according to the simple formula given above. This value is close to 
£^F,crit given in Eq. fl36l) and thus one could expect a strongly reduced repulsion or even 
a net attraction between pairs of these glass spheres. A corresponding extension of this 
single-sphere experiment would be highly desirable. 

For truly nanoscopic colloids {tq < 1 fim) with this same charge density on the air or 
oil side, the formula above predicts ep = 0(10^^) and therefore the electrostatic repulsion 
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would always dominate the capillary attraction. This is consistent with the results in Ref. 
obtained for polystyrene spheres {tq ^ 1 /im) at the oil/water interface. For highly salty 
water, the charge density on the colloid/oil interface is estimated experimentally to be 
(T ~ 20 /iC/m^ and force measurements between two spheres confirmed the repulsive dipole- 
dipole interaction with no sign of capillary attraction. 

Assumption (B) amounts to modelling a system in which the Debye length of water 
is comparable or larger than the radius tq. We have studied in detail the case that the 
dielectric constant of one fluid phase (e.g., water) is much larger than the dielectric constant 
of the other fluid phase (e.g., air). The capillary-induced potential V^en is asymptotically 



This work is de facto a correction 24] of previous work 19|. 
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attractive (see Fig. Hj) but the total potential Vmen + K-ep can be asymptotically attractive 
for values £i? ~ 1 only if /tro > 1 (see Fig. [5]). However, unlike the model corresponding 
to assumption (A), even in such a case Vmen + Kep becomes repulsive at small separations 
(Fig. [H]). This effect can be traced back to a crossover in the interfacial stress n(r) from the 
asymptotic algebraic decay (as in model (A)) to an exponential decay at closer distances 
r < 7/k from the particle (see Fig. [3]). 

In contrast to the experiments carried out with nano- and microcolloids at interfaces with 
salty water, some experiments j^, [lo|, 12, [isl, Q] have been performed with microcolloids 



at interfaces of ultrapure water such that ~ ro ~ 1 /im. These experiments could 
have explored phenomena beyond the small deformation regime. Equation (1471) provides a 
relationship between ep and the relevant parameters of the experimental system; regretably, 
the value of the total charge q is usually uncertain. In terms of the surface charge density 
(J one has q = 27ro"(ro sin 6')^(1 + cos^), because for ultrapure water the electrostatic field 
is dominated by the unscreened charge of the particle on the water side. For typical values 
€2 = 81, 7 = 0.05 N/m, n = 1 /im~^, and 6 = 7t/2 (so that tq = radius of the particle), 
Eq. ( H7|) gives ep ~ 371 cr'^ r^Vlro) with tq in fim and a in units of e/nm^. The values for a 
quoted in the literature range from 0.07e/nm^ (Ref. Q) to 0.53 e/nm^ (Ref. [s]). Accordingly, 
Fig. ([7j) shows that it seems possible to have capillary attraction (i.e., ep > ^F.crit ~ 1) for 
typical values of the particle radius in the micrometer range. However, we emphasize that, 
unlike the conclusion concerning the asymptotic decay of Vmen, the expressions relating the 
value of Ep and ep^cnt with the parameters of the system involve the behavior of the electric 
field near the particle (see Eqs. (H8ll and (l53l) ). Therefore, they are expected to be affected 
by corrections to the monopolar approximation and the Debye-Hiickel approximation 35|. 
Thus, further theoretical work is required to understand these experiments properly and to 
arrive at reliable predictions. 
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FIG. 7: Parameter space spanned by the surface charge density a (m units of e/nm^) and the 
contact hne radius tq (in ^m) (see text for the fixed values of the other parameters). The solid line 
corresponds to the loci ep = eF.criti so that capillary attraction is predicted to occur in systems 
the parameters of which fall into the region above this curve. As a reference curve, the dashed line 
corresponds to the loci ep = 1. 




APPENDIX A: INTERFACE DEFORMATION FIELDS 

In this appendix we derive the single-colloid deformation u{r) and the correction ^^(r) 
to the superposition approximation for two colloids. 

The deformation field u{r) in the single-colloid configuration is readily obtained as the 
rotationally symmetric solution of Eq. ([5]): 

L 1 s 

u(r) = ro(en — £f)^^ / (issn(s)ln- 

r 1 Jr r 

r^-" if £p-£^ = 0, 

(ro < r < L and n > 2). (Al) 

Inr if Ep — Eu ^ 0, 



The limit L — > oo is well defined if Ep — Eu = 0] otherwise, the presence o 



the boundary 



condition (ISc|l is required to regularize the possible logarithmic divergence [2^ . The physical 
interpretation of this regularization is a force acting on the boundary Cl of the interface 
which compensates the net force 27i'yrQ{Ep — En) localized around the particle. Since these 
two forces act at well separated locations, there is an intermediate range of lengths where 
there is approximately no force acting on the interface, so that the corresponding deformation 
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FIG. 8: The position of a point P near a colloid is parametrized by the polar coordinates (p, 97). 
The unit vector normal to the (circular) contact line is n = ep. 



varies logarithmically. (The electrostatic analogy developed in Ref . l27| provides a transparent 
visualization of this explanation.) 

The correction Mm(r) to the superposition approximation cannot be computed analytically 
in an easy manner. But we shall derive several asymptotic properties of the solution as the 
interparticle separation d becomes large. In order to follow the arguments, the reader will 
find the electrostatic analogy useful, which relies on the formal analogy of Eq. ([8]) with the 
equations for the two-dimensional electrostatic potential: Um plays the role of the potential 
and 2Ilm the role of the charge density. (This analogy is worked out in detail in Ref. 1271 .) 

First we approximately compute the function Um{r) near the colloids, i.e., at distances 
fQ^P'-= |r— rQ,| <^ d from the particle (see Fig. [8]). Concerning the boundary condition fl8bp 
at p = To, we note that 



-rin 



Up - (up) 



£12 - 



ro 



nn : VVu{d) - -V^u{d) 



u"{d) 



u'{d) 
d 



COs2(y9, 



(A2) 



which follows from expanding the single-particle solution (lAip around r = d, where one 
finds (see Eq. flAT]) ) 

u'{d) [d-"", if£^-£n = 0, 
d-\ 



u"{d) 



(A3) 



d I Hep -en 7^0. 

The field Umir), according to the model discussed at the beginning of Subsec. f ill Bp . is 



peaked at the colloids and decays far from them, so that it has approximately rotational 
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symmetry in the range tq < p d. More precisely, since 11^ is proportional to |n(|pep + 
dex\)Ii{p)\^^'^, one finds, after expanding in terms oi p/d <^ 1 and in view of Eq. f|T3|) . 

n^(r) ~ |n(d)n(p)|V2 ^ 0(^-n/2-i)_ ^^4^ 

Accordingly, in leading order in l/d Eq. ([8]) turns into 



pdp 
n ■ Vu„ 



dum 



dp 



+ ^ --^M - d--'\ roKp-^d, (A5a) 

p^ o^p^ 7 



£^n„ ~ d if - en = 0, 

jp=ro I -(ro/4)[u"(d) /(i] cos 2(y9 ~ ci"^ if e^r - 7^ 0. 

(A5b) 



Thus, if — en = the near-particle solution is approximately rotationally symmetric and 
dominated by the pressure field 11^: 



if £_p — £n = 

' ro P 



p 2 fP 

Um{p) ~ Ao + ro£:n,„ln h- l ds sYim{s)\n.- 

ro 7 Jro 

2 f°° s 2 f°° p 

^ Aq-I — / ds sllm{s) In 1 — / ds s Ilm{s) In - , (A6a) 

7 Jro ro -fjp s 

where the second expression follows by inserting the estimate (|T7j) for £n,„- Only the last term 
depends on p, while the second one is an additive constant proportional to ^yil{d) ~ rf^"/^. 
On the other hand, if ei? — en 7^ the near-particle solution is dominated by the single- 
particle solution for the boundary condition: 



Umip) ~ Ao+- 



, , u"(d) + -A2 
12 \ d ^ W 3 



2 / \ 2 



rl ( u\d) \ 1 ^ ( rQ\ , , ( p 



+ — > cos 2^9, if e^-en 7^ 0. 



(A6b) 

The integration constants Ao and A2 are determined by the solution far from the particles 
and for our purposes here, it suffices to provide an estimate of how they depend on the 
separation d. In view of the electrostatic analogy the solution Um far from the particles, 
r ^ d, can be expressed in terms of a multipolar expansion. The "capillary monopole" 



Q^:= dA2Ii^ + -i <h din-{- Vu^) (A7) 



is found to vanish exactly by virtue of the boundary condition ( ]8bl) . The "capillary dipole" 
vanishes due to the refiection symmetries of the configuration shown in Fig.[2]upon X — > —X 
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and y — i> —Y. But in general the "capillary quadrupole" will be nonzero. Therefore 
the distant field is (C is a proportionality constant) 

d r 

Um{r>d)r^^ + C\n- / dssUmis), (A8) 

^ ^0 Jr 

where we have neglected the angular dependence of the quadrupolar field (we are interested 
only in the decay with distance) and dm denotes the typical value of the elements of the 
quadrupole Dm- The second term is approximately the "potential" created by the "capillary 
charge" beyond r, which is not accounted for by the multipolar moments, and arises as a 
correction to the multipolar expansion because the "capillary charge" 11^ does not have a 
compact support. It scales as r^"*^ Inr. On the other hand, because the "capillary monopole" 
of a single particle is of the order of ■jroeu^, the quadrupole of the "capillary charge" 
distributed over a region of size ~ d will be dm ~ d'^'yvQEn^ ~ d'^ Thus, the second 
term is negligible compared to the quadrupolar term if r ~ rf, provided n > 4, so that one 
finally finds Umir ~ c?) ~ For reasons of consistency we expect that the near-particle 

solutions in Eq. flA6l) should scale like this if extrapolated to p ~ d: if Eq. (IA6al) is evaluated 
at p ~ c/ one obtains Aq = 0(rf""/^), while from Eq. (lA6b|) it follows that A2 = 0{d-^). 

In sum, the amplitude of the near-particle solution scales with the interparticle separation 
as follows: 

frf-"/^ ifeF-£n = 0, 
Mm(p~ro)~< _ (n>4). (A9) 

[ d-\ if ep-euT^ 0, 

APPENDIX B: FLOTATION FORCE AND DISJOINING PRESSURE 

In this appendix we discuss briefly how our results are modified if the gravitational force 
is relevant. We also show how the same formal results hold if the lower fluid phase is a thin 
film on which, instead of gravity, dispersion forces due to a confining substrate are acting, 
exerting the so-called disjoining pressure. 

1. Flotation force 

The effect of the acceleration of gravity g gives rise to an additional contribution to the 
free energy ([2]): the gravitational potential energy of the fluids with respect to the reference 
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FIG. 9: Schematic drawing of a fluid film on top of a substrate. I is the thickness of the film in 
the reference, flat configuration, and u is the deformation of the interface. 



configuration is 



1 r / - \ 2 
1 , . I u 



where we have introduced the capillary length A := A/7/[(f?- — Q+)g] in terms of the mass 
densities g+ and of the upper and lower fluid phases, respectively. This length has typical 
values in the millimeter range. A pressure field (force per unit area) can be associated with 
this free energy: 



n.^-W- -i^= -i^W- (B2) 



6^ a;rav T 

5u{r) ~ ~ A2 

This expression can be inserted directly into Eq. (I3b|) . and one finds that at large distances 
from the particles u{r) ~ exp(— r/A) and the field Hgravlr) does indeed decay sufficiently fast. 
However, one cannot apply the results we have derived previously without certain changes 
because Ilgrav depends explicitly on u. Nevertheless, in the limit of a large capillary length 
one can neglect J?-grav (and thus Ilgrav) altogether and retain only the gravitational force 
acting directly on the colloidal particle, i.e., en = and F = F is the weight of the colloidal 
particle (corrected for buoyancy effects). In this case Vmen{d) is given by the superposition 
approximation (see Eq. (12T|) ) and reproduces the flotation force in the regime d -C A jisl . 



2. Disjoining pressure 

If the size of the particle lies below the micrometer, gravity is quantitatively negligible 



22] . The same formalism, however, is applicable in the experimentally relevant case that the 



lower fluid phase is a film of thickness h{r) = i + u{r) on top of a (solid or liquid) substrate 
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(see Fig. [9]). If £ is within the range of the underlying dispersion forces, an additional 



contribution JFjigp to the free energy ([2]) arises 
term) : 



(neglecting a constant, u-independent 



^ Hi 



disp 



dA 



H 



+ (£ + M)AnA/i 



(B3) 



where H is known as the Hamaker constant. An is the number density difference between 
the bulk phases the film and the upper phase belong to, and A/x is the undersaturation of the 
upper phase in terms of the chemical potential. The condition that the flat film of thickness £ 
is an equilibrium solution imposes the relation An A/x = IJijP' . The corresponding pressure 
field associated with JF^jsp is called "disjoining pressure": 



Ildisp(r) 



(B4) 



In the regime of small deformations one has |m| ^ £ and J?-disp can be expanded around 
M = yielding 



^disp ~ ^disp[M = 0] + -7 



dA 



u 



A, 



disp 



(B5) 



which has the same form as Eq. ( IBll) with an effective "capillary length" Ad isp = Vt^V(6^7 
(which actually coincides with the so-called lateral correlation length 36]). Thus Hdisp 
takes the form of Eq. flB2p with A replaced by Adisp- H the film is sufficiently thin, the 
length Adisp can be so small that it becomes relevant and a phenomenology may arise which 
is similar to the macroscopic one induced by gravity. For typical values H = 10^^° J and 
7 = 0.05 N/m, one has Adisp = 1 /^m for a film thickness i = 0.033 /xm. 



APPENDIX C: EFFECTIVE CAPILLARY POTENTIAL ENERGY IN THE 
CASE OF AN IDEALLY CONDUCTING FLUID 

In this appendix we provide the mathematical steps leading to Eq. (!32|) which is valid in 
the limit d ^ and from which the effective capillary potential energy is obtained. 

In order to compute en™ , "we insert the ansatz (!3T|) into the definition (]15bp . In the limit 
d oo, the main contribution to the integral stems from the regions around each particle 
so that 

2 

£n™~— v/n(rf)/ rfrr^nW. (CI) 
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With the pressure field given by Eq. (129|) this reduces to 



(C2) 



and Eq. (l32k ) is obtained upon performing the integral (Eq. 3.191.2 in Ref. |32| ). 

For the evaluation of the integral in Eq. fl32b ) we use the solutions flAll) and flA6aP valid 
ioT ep- eu^ = 0. With Eq. ([29]) this leads to 



Vm = — ds s n(s) = erb(fi)eF — / c/x (x — 1 

J r Jr/ro 

= ^erb{fi)eF (^j) 2i^i (l - 4; 5; ^) 



(C3) 



and, with the ansatz (I3T1) . to 



7r 



Jr T \ d / Jr/ro 



2 r 



(C4) 



(Concerning the last lines in Eqs. (1C3P and (lC4p see Eq. 3.194.2 in Ref. 
one to obtain 



32.) This enables 



I dA (Vm„) ■ (Vm2) ^ nh{iifel (^)' x 

'S'mcn,2 



(C5) 



X I dr r 



rt3 



To 



2F1 l-/i,4;5;^ 2F, 



2 r 



so that Eqs. fl32b . 1551) follow upon a change of variable in the integral. 



APPENDIX D: ASYMPTOTIC BEHAVIOR OF THE ELECTROSTATIC PO- 
TENTIAL AND OF THE ELECTRIC FIELD 

In this appendix we discuss the asymptotic behaviors of the functions 1n{k) defined in 
Eq. fj43l) . First we mention that at first sight the integrals may appear to be divergent due 
to the weak power law decay of the integrands, which are, however, oscillatory. Actually, 
the integrals are regularized by an exponential, so that for instance 

/"°° X 
Uk):=]im dxe~''^Mx)^^—^===, (Dl) 
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reflecting the physical situation that the charge q is positioned at a (dimensionless) height 
h > above the flat interface. In the following mathematical manipulations this regulariza- 
tion scheme is implied, which, unless required, we do not write explicitly to avoid a clumsy 
notation. 

I — I 

Reference provides a method to obtain the expansion of 2a{k) in powers of ei/e2 <^ 1. 
The idea is to split the integrals into a sum of two terms involving only odd or even powers 
of the ratio ei/e2, respectively. This leads to^ 



J„(A;) = e"'^ + ^ -k{Io{k)-Lo{k))-l + ^(ei/es)', k<^{e2/e^)^\ (D2) 

62 L/ -I 

in terms of the Bessel function Jq and the Struve function Lq . We note that the validity 
of this expression is restricted to sufficiently small values of k. However, for the typical values 
of the ratio 62/61 occurring in the experiments so far, this does not impose any physically 
relevant constraint on k. The asymptotic behavior for large k is 

J,(fc)^e-'= + ^-3-, 1 « « (62/61)2. (D3) 
62 k^ 

As a general result the coefficients of the odd powers of 61/62 decay algebraically for 
k ^ 1, while the coefficients of the even powers decay exponentially. With Eq. (H3b ) this 
leads to 

X,(A;) ^ A;e-^' + 3--^^, 1< A; < (62/61)^ (D4) 
62 k^ 

The procedure employed in Ref.y can be extended to analyze Itik). The radical is eliminated 
from the denominator of the integrand in Eq. (H3b) and the integral is split as follows: 



62 



00 



^b{k) = — dx Jo X / N21 2 I /2 ~ / M / / N21 2 I /2 - 

61 Jo [1 - {(^i/(^2r\x^ + k^ Jo [1 - {(^i/(^2>^\x^ + k^ 



The first term involves only odd powers of 61/62, the second term only even powers. The 
first integral can be rewritten as 



62 r, x{x' + e) 62/61 r . 

— / dx Jo X / N21 2 I /2 = 1 rn~^ / ^ 



k 



2 fe2 



^ [l-(eiA2)^] 



[l-(^lA2)2] 

(D6) 



^ We note (i) a misprint in Eq. (10) in Ref. 0, because there the term 2™ should be in the numerator, and 
(ii) that the series provided in that Eq. (10) can actually be carried out and is proportional to the error 
function. 
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With the identities (see Eqs. 6.623.2 and 6.532.4 in Ref. M) 



and 



dx Jn(x) X = 



X 

dx Jo{x) ^ ^ = Ko{a) 



the first integral in Eq. ( ]D5I) can be written as 



€2 J, Xjx^ + k') 

ei Jo '^"'^ [1 - {e,/e,y]x^ + F 



ei/e2 



[1 - (ei/62; 



212 



(D7) 



(D8) 



1 - ieje.y) 



The second integral in Eq. fIDSP can be written to leading order in ei/e2 as 



X 



2 + k^ 



dxJoix) 
+ fc^ 



dx Jnix) 



fc2 



(D9) 



(DIO) 



We introduce the exponential regularization as exp {—hyx^^+W) and note the identity 
(Eq. 6.637.1 in Ref . Id) 

'k 



00 ^-h\/x'^+k'^ 

dx Jq{x) = 
yx^ + A;^ 



Therefore, the righthand side of Eq. (1D10I) can be rewritten as 



dx Jq[x) „ , ,„ — = lim 

x^ + k^ h^O 



dx Jo{x) = 
vx^ + A;^ 



(DU) 



(D12) 



Therefore from Eqs. flD9p and flD12p one finally arrives at the following expansion: 
1 



Uk) = 



^ e Ko{k) + 0{e,/e,f. (D13) 
£2 



This expansion may also be restricted to sufficiently small values of fc, although we have 
not found this upper bound as function of ei/e2. However, as argued above this possible 
constraint is expected to be physically irrelevant because typically ei/e2 is sufficiently small. 
The asymptotic behavior of this expression for ^ 1 is 



1 < k, e2/ei. 



(D14) 
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